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OBJECTIVES 

There  has  been  significant  interest  in  applying  geometric  techniques  from  dynamical  systems  to  the 
study  of  transport  and  mixing  in  oceanic  flows.  Much  of  this  has  been  based  on  the  idea  that  stable  and 
unstable  manifolds  of  hyperbolic  trajectories  provide  a  useful  geometric  description  of  transport  and 
mixing  processes  when  the  velocity  fields  are  two-dimensional  and  periodic  in  time.  Original  objectives 
of  this  research  effort  included: 

•  applying  these  dynamical  systems  techniques  to  flows  with  aperiodic  time-dependence 

•  identifying  additional  diagnostics  relevant  to  transport  that  can  be  acquired  fi-om  the  invariant 
manifold  template 

•  extending  the  functionality  of  the  existing  yftool  numerical  package  and  improving  the 
computational  efficiency  and  numerical  accuracy  of  the  algorithms 

•  making  the  computational  tools  more  accessible  to  others  including  better  interfaces  for 
incorporating  new  models  and  integration  with  computational  environments  such  as  Matlab 

The  software  is  especially  suited  to  analyze  velocity  fields  that  are  given  as  data  sets;  velocity  fields 
derived  either  as  output  from  numerical  simulations  of  the  governing  dynamical  equations,  or  derived 
from  observations  in  the  real  environment. 

BACKGROUND 

Within  time-dependent,  two-dimensional  velocity  fields,  regions  of  strong  saddle-type  dynamics  play  a 
significant  role  in  enhancing  transport  and  mixing  of  fluid.  When  the  time-dependence  is  periodic, 
asymptotic  {t  ->  ±oo)  structures  such  as  stable  and  unstable  manifolds  of  hyperbolic  periodic 
trajectories  provide  a  geometric  template  describing  this  transport  and  eventual  mixing.  In  this  case  it  is 
sufficient  to  analyze  the  dynamics  of  the  two-dimensional  Poincare  map,  where  the  transport  is 
described  by  the  transverse  intersections  of  one-dimensional  stable  and  unstable  manifolds  associated 
with  hyperbolic  fixed  points  [5],  [9]. 
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When  the  velocity  fields  exist  as  data  sets,  the  time-dependence  is  generally  aperiodic,  the  fields  are 
defined  only  on  a  finite  interval  in  time,  and  a  single  two-dimensional  map  is  no  longer  sufficient  to 
describe  the  transport  dynamics.  This  aperiodic  time-dependence  and  finite-time  aspect  of  numerical 
velocity  fields  motivated  new  computational  techniques  for  identifying  dynamical  structures  useful  in 
characterizing  transport  and  mixing.  Methods  developed  with  collaborators  at  Brown  University  and 
Woods  Hole  are  used  to  define  certain  distinguished  two-dimensional  invariant  manifolds  which  behave 
much  like  asymptotic  stable  and  unstable  manifolds  on  the  available  finite  time  interval  [7],  [3].  These 
effective  invariant  manifolds  can  be  identified  numerically  by  evolving  an  appropriate  set  of  passive 
tracers  (solutions  to  the  ODE,  x'(0  =  v(t ,  x)).  If  the  initial  tracer  positions  straddle  the  hyperbolic 
region,  the  fluid  particles  trace  the  unstable  manifold  in  forward  time  and  the  stable  manifold  when  the 
ODE  is  run  in  reverse  with  respect  to  time.  As  the  manifolds  lengthen,  the  numerical  algorithm  must 
start  additional  tracers  to  keep  the  manifolds  accurately  resolved. 

When  viewed  at  a  fixed  time,  the  manifolds  for  two-dimensional  vector  fields  appear  as  curves  in  the 
phase  space.  Segments  of  these  manifolds  can  be  used  to  define  time-dependent  boundaries  (Lagrangian 
boundaries)  between  different  regions  of  the  flow.  When  the  manifolds  intersect,  lobes  are  formed; 
closed  regions  in  the  phase  space  bounded  by  two  segments  of  manifolds,  one  stable  and  the  other 
unstable,  connecting  a  pair  of  intersection  points  (hetero-  or  homoclinic  trajectories).  As  these  lobes 
evolve  in  time,  they  are  associated  with  an  exchange  of  fluid  between  the  two  regions,  a  mechanism 
often  referred  to  as  “chaotic  advection”  or  “chaotic  transport”.  The  area  of  the  lobes  formed  from  the 
intersection  of  these  manifolds  can  be  used  to  quantify  the  material  transport  across  these  boimdaries  [7]. 

RESULTS 

Barotropic  meandering  jets 

In  work  with  Rogerson,  Pratt,  and  Jones  [10],  stable  and  unstable  manifolds  were  identified  in 
simulations  on  a  fully  nonlinear,  dynamically  consistent,  barotropic  model  of  a  meandering  jet.  The 
manifolds  confirmed  that  the  most  significant  exchange  of  fluid  is  between  the  eddies  that  bound  the  jet 
and  the  region  outside  the  jet.  Fluid  exchange  between  the  eddies  and  the  jet  itself  was  4-7  times  smaller, 
consistent  with  earlier  results  involving  simpler  kinematic  models  [6].  When  dimensionalized  to  the 
scale  of  the  Gulf  Stream,  the  voliune  transport  associated  with  this  mechanism  of  "chaotic  advection" 
was  comparable  to  the  fluid  crossing  the  Gulf  Stream  due  to  the  mechanism  of  ring  detachment. 

In  a  related  result  by  Lozier  et  al  [4],  RAFOS  float  trajectories  in  the  vicinity  of  the  Gulf  Stream  were 
viewed  in  a  moving  frame.  The  fate  of  the  trajectories  were  then  interpreted  based  on  the  dynamical 
systems  analysis  from  this  numerical  model  of  meandering  jets. 

Viscous  perturbations  of  vorticity  conserving  flows  and  separatrix  splitting 

In  joint  work  with  Sandstede,  Balasuriya,  and  Jones  [II],  these  numerical  techniques  were  used  to 
identify  stable  and  unstable  manifolds  in  a  time-dependent,  two-dimensional  velocity  field  that  satisfies 
a  viscous  perturbation  of  the  barotropic  vorticity  equation  on  the  )8-plane.  This  paper  extended  earlier 
results  from  Balarsuriya’s  Ph.D.  thesis  [1]  by  developing  a  Melnikov  theory  applicable  to  vector  fields 
defined  over  finite,  but  sufficiently  large  time  intervals.  In  this  model  problem  we  had  an  exact  solution 
to  the  perturbed  barotropic  equation.  The  numerical  computations  confirmed  the  predicted  separatrix 
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splitting  and  the  non-existence  of  heteroclinic  points  (see  Figure  1).  The  distance  of  the  splitting  as 
determined  from  the  numerical  simulations  compared  favorably  with  the  theoretical  values. 


Exact  Solution  to  Barotropic  PV  Equation  (e  =  1 .0-10"^ ) 
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Figure  1.  Computed  manifolds  for  e  =  1  x  lO"!  The  particle  trajectories  satisfy  the 
equations,  x'it)  =  — y,  i)  and  y(/)  =  Pi  t),  where  the  streamfunction 

ijr  satisfies  the  partial  differential  equation,  qt  -  ^yqx  -I-  "ifxqy  =  eA^  with  potential 
vorticity  q  =  A\|f  +  ^y. 


Chaotic  transport  for  an  island  recirculation 

In  work  with  Pratt,  Jones,  and  Helfnch  [12],  lobe  analysis  was  used  to  analyze  a  time-dependent, 
boundary-trapped  recirculation.  The  recirculation  gyre  occurs  in  a  numerical  model  of  wind-driven  flow 
around  an  island,  with  a  geometry  similar  to  persistent  eddies  such  as  the  Alboran  Gyre  and  The  Great 
Whirl.  In  the  numerical  model,  an  anti-cyclonic  circulation  is  driven  by  a  steady  wind  stress,  setting  up 
a  large  recirculation  gyre  to  the  east  of  the  island.  The  forcing  is  parametrized  as  the  ratio  of  the  inertial 
boundary  layer  thickness,  8;  =  (t///3) to  the  Munk  boundary  layer  thickness,  Sm  = 

Here  U  denotes  a  horizontal  velocity  scale  for  the  flow  and  A  h  the  horizontal  eddy  viscosity. 
Throughout  the  numerical  experiments  Sm  was  held  fixed  at  40  km.  For  =  1 .0  the  flow  is  steady 
and  the  only  material  transport  between  the  gyre  and  the  surrounding  fluid  is  due  to  a  downward 
“Ekman  pumping”,  a  volume  flux  that  can  be  estimated  from  theory. 

As  Sj/Sm  increases,  the  two-dimensional  velocity  field  becomes  time-dependent  although  the 
recirculation  attached  to  the  island  boundary  persists  for  long  times  providing  an  ideal  coherent 
structure  for  studying  Lagrangian  transport.  The  effective  stable  and  unstable  manifolds  emanate  from 
hyperbolic  regions  close  to  the  island  boundary  and  at  either  end  of  the  gyre.  As  the  forcing  is  increased 
in  these  simulations,  the  lobe  geometry  becomes  less  regular  in  both  time  and  space.  Representative 
lobes  are  shovm  in  Figure  2  for  increasing  values  of  Sj/Sm- 
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Figure  2.  Lagrangian  transport  for  the  numerical  simulations  of  flow  past  an  island. 
The  Munk  layer  thickness,  bm,  is  fixed  at  40  km  in  all  four  runs.  The  manifolds  were 
evolved  for  50  days.  Te  denotes  the  Ekman  transport  and  Tc  the  chaotic  transport. 


The  analysis  shows  that  the  transport  out  of  the  gyre  due  to  Ekman  pumping  {Te)  is  quickly  surpassed 
by  transport  fluxes  associated  with  chaotic  advection  (Tc).  Figure  3  is  a  graph  comparing  the  Elman 
transport  to  the  transport  due  to  chaotic  advection.  The  Ekman  transport  is  equal  to  the  downwards 
Ekman  pumping  velocity  integrated  over  the  fixed  area  of  the  recirculation  as  determined  by  linear 
theory  [8].  A  least-squares  fit  yields  the  parametrization,  Tc  *=»  0.54[((5//5;i/)^  —  (1.8)^]^-^.  The  shaded 
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region  in  Figure  3  is  an  estimate  of  the  range  xtiSj /8m  over  which  the  method  of  lobe  dynamics  and  the 
Lagrangian  definition  of  the  recirculation  boundary  is  applicable  and  helpful. 


transport  (Sv.) 


subjective  limit 
of  utility 


(8i/5m)2  (=2.752) 

Figure  3.  Ekman  transport  (solid  line)  vs.  transport  due  to  chaotic  advection  (dashed 
line).  Note  that  Te  is  independent  ofhM  and  is  regarded  here  as  a  function  of&r/^M 
only  insofar  as  is  fixed  in  the  numerical  experiments. 

Potential  vorticity  budgets 

As  part  of  the  analysis  for  the  island  recirculation  we  investigated  methods  for  characterizing  the 
potential  vorticity  budget  of  the  recirculation  gyre.  We  were  interested  in  formulating  approaches  that 
consider  the  Lagrangian  definition  of  the  gyre  as  compared  with  a  fixed  recirculation  boundary  defined 
from  the  time-averaged  velocity  field.  The  most  promising  approach  defines  the  Lagrangian  boundary 
as  segments  of  stable  and  unstable  manifolds  joined  by  a  "gate".  With  this  representation  of  the  gyre, 
mass  transport  takes  place  through  the  gate  continuously  in  time.  The  analysis  showed  that  the  fixed 
boundary  firom  the  time-averaged  flow  exaggerates  the  contribution  fi-om  advection  of  vorticity  into  the 
gyre.  This  appears  to  be  due  in  part  to  the  recirculation  passing  as  a  whole  across  the  fixed  boundary, 
rather  than  a  true  vorticity  flux  at  the  boundary  of  the  gyre. 
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Advances  in  software  (vftool) 

This  effort  included  a  number  of  significant  improvements  to  the  numerical  tools  used  for  performing 
the  dynamical  systems  analysis. 

•  added  adaptive-step  method  for  solving  the  differential  equations  (Runge-Kutta-Fehlberg  method) 

•  overhauled  package  of  C-code  to  simplify  the  model  inputs,  improve  the  command-line  interface, 
and  add  help  files 

•  code  now  runs  under  irix,  linux,  netbsd  operating  systems 

•  overhauled  Matlab  toolboxes  for  performing  data  analysis,  eliminated  i/o  problems  due  to 
system-dependent  byte  order  for  unformatted  data,  added  complete  help  pages  for  all  m-files 

Several  other  programming  features  also  received  attention  but  had  not  reached  the  point  of  being  fully 
implemented  into  the  working  code  by  the  end  of  the  funding  period. 

•  solving  the  linearized  variational  equations  along  with  the  nonlinear  differential  equations;  this 
will  provide  additional  characterization  along  the  lines  of  Haller’s  work  [2] 

•  parallelizing  the  code  to  speed  up  manifold  calculations  which  require  tracking  thousands  of  fluid 
particles 

•  representing  the  material  curves  using  a  C*  interpolation;  this  includes  tracking  derivative 
information  using  the  linearized  flow  and  using  error  estimates  to  represent  the  curves  more 
efficiently 


PERSONNEL 

In  addition  to  the  P.I.,  this  program  was  supported  by  the  efforts  of  one  graduate  student,  Paul  von 
Dohlen,  and  several  undergraduate  students,  William  Beksi,  Shilpa  Savaliya,  Sunil  Shah,  and  Wei  Weng. 


6 


REFERENCES 


[1]  H.S.  Balasuriya,  C.K.R.T.  Jones,  and  B.  Sandstede,  1998:  Viscous  perturbations  ofvorticity 
conserving  flows  and  separatrix  splitting,  Nonlinearity  1 1, 47-77. 

[2]  G.  Haller,  2000:  Finding  finite-time  invariant  manifolds  in  two-dimensional  velocity  fields,  Chaos, 
10(1):99-108. 

[3]  G.  Haller  and  A.  Poje,  1998:  Finite-time  transport  in  aperiodic  flows,  Physica  D,  1 19(3-4):352-380. 

[4]  M.S.  Lozier,  L.J.  Pratt,  A.M.  Rogerson,  and  P.D.  Miller,  1997:  Exchange  geometry  revealed  by  float 
trajectories  in  the  Gulf  Stream,  J.  Phys.  Oceanogn,  27, 2327-2341. 

[5]  R.S.  MacKay,  J.D.  Meiss,  and  I.C.  Percival,  1984:  Transport  in  Hamiltonian  systems,  Physica  D, 
13:55-81. 

[6]  P.  Miller,  C.  Jones,  G.  Haller,  and  L.  Pratt,  1996:  Chaotic  mixing  across  oceanic  jets.  Chaotic, 
fractal,  and  nonlinear  signal  processing,  AIP  conference  proceedings  375,  R.  Katz,  editor,  591-604. 

[7]  P.D.  Miller,  C.K.R.T.  Jones,  A.M.  Rogerson,  and  L.J.  Pratt,  1997:  Quantifying  transport  in 
numerically-generated  velocity  fields,  P/iy5fca  D,  110, 105-122. 

[8]  J.  Pedlosky,  L.J.  Pratt,  M.A.  Spall,  and  K.R.  Helfrich,  1997:  Circulation  around  islands  and  ridges, 

J.  Marine  Res.,  55, 1 199-1251. 

[9]  V.  Rom-Kedar  and  S.  Wiggins,  1990:  Transport  in  two-dimensional  maps,  J  Arch.  Rat.  Meek 
Anal,  109,  239-298. 

PUBLICATIONS 

[10]  A.  Rogerson,  P.  Miller,  L.  Pratt,  C.  Jones,  1999:  Lagrangian  motion  and  fluid  exchange  in  a 
barotropic  meandering  jet,  J.  Phys.  Oceanogr.,  28, 2635-2655. 

[1 1]  B.  Sandstede,  S.  Balasuriya,  C.  Jones,  and  P.  Miller,  2000:  Melnikov  theory  for  finite-time  vector 
^Q\ds,  Nonlinearity,  yo\.  13, 1357-1377. 

[12]  P.  Miller,  L.  Pratt,  K.  Helfrich,  and  C.  Jones,  2002:  Chaotic  transport  of  mass  and  potential 
vorticity  for  an  island  recirculation.  Journal  of  Physical  Oceanography,  vol.  32,  no.  1,  80-102. 


7 


